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Abstract 

A systematic mathematical methodology for derivation of boundary layer expansions 
is presented. An explicit calculation of boundary layer sizes is given and proved to be co- 
ordinates system independent. It relies on asymptotic properties of symbols of operators. 
Several examples, including the quasigeostrophic model, are discussed. 



1 Introduction and main result 

Boundary layers appear in various physical contexts, such as fiuid mechanics, thermodynam- 
ics, or ferromagnetic media. From the mathematical point of view, they are related to singular 
perturbation problems in bounded domains. The singular perturbation is due to the presence 
of small parameters in the dimensionless governing equations. For instance, in magnetohy- 
drodynamics (MHD), the so-called Rossby number (ratio between the angular velocity of the 
fiuid and the angular velocity of the Earth) can be as low as 10~^. Similarly, the Prandtl 
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number (ratio between hydrodynamic and magnetic diffusivity) is about 10~^ in the Earth's 
core. In the interior of the domain, these small parameters lead to some reduced dynamics of 
the equations. For instance, in highly rotating fluids, the velocity field does not vary along the 
rotation axis, see textbook j^. This reduced dynamics is often incompatible with boundary 
conditions. This yields boundary layers, in which the solutions of the equations have strong 
gradients, in order to satisfy the boundary conditions (typically a Dirichlet condition at a 
rigid surface). 

In a formal setting, boundary layer problems are connected to systems of the type 

A'{x,D,X + Q"{u") = f, xen, (1.1) 

where e G R"' describes the small parameters of the system, Q C is the domain, u^(x) 
the unknown, /^(x) the data, A"^ is a linear differential operator (most often of elliptic or 
parabolic type), and the nonlinear part. One must of course supply these equations with 
appropriate boundary conditions. 

The basic idea, which underlies the mathematical study of all boundary layers, is that the 
solution of (|1.1|1 should satisfy an asymptotics of the form 

u%x) ~ .(X) + < (x, ^^^^) + ... + uf (x, ^^^^) , e^O, 



ai(e) J V ' ar(e) / ' ' (1.2) 

(^i) ^0, 6*^^+00, Vi 

This means that should have a regular part, depending on x, but also a singular part, 
depending on stretched variables 9i = d{x, dQ,)/ai{e), with ai{e) — > as e ^ 0. This singular 
part should be localized near the boundary and express the strong gradients of boundary 
layers. Broadly speaking, the aim of boundary layers studies is to answer the following 
questions: 

i) What are the possible ai{e) (the boundary layer sizes)? 

ii) What are the possible profiles u, uf ? 

iii) Is the asymptotics (|1.2|) correct (in a sense to be determined) ? 

The first two questions are related to the derivation of boundary layers, whereas the third 
one is connected to stability issues. 

The aim of this note is to give some insight into the derivation problem. In some cases, the 
derivation is quite easy, and the stability analysis is the most difficult part (see for instance 
|1U1 115j on viscous perturbations of hyperbolic systems) . But in most situations of physical 
interest, it may involve a variety of length scales and equations. A typical example is the 
description of water in a highly rotating tank. If the water is at rest in the rotating frame 
attached to the tank, it is well modeled by Stokes equations with Coriolis term: 

e X u + Vp — EAu = f, , , 

(1.3) 

V-u = 0, 

where e = (0, 0, 1) is the rotation vector, and E is a small parameter called the Ekman 
number. For such system (with appropriate source term / and boundary conditions), the 
structure of the solutions near a flat horizontal boundary is understood: there is a boundary 
layer of size E^^'^, the Ekman layer (see [Oj). In this case, one can even carry an analysis 
of the full Navier-Stokes equations with rotation: we refer to jl21 1141 0] among others. But 
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Figure 1: Boundary layers of rotating fluids, near a sphere and at the circumscribing cyUnder 
(fohowing Stewartson ^ 
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in more elaborate geometries, following the articles of Stewartson I19j. many other layers 
develop: for instance, between two concentric spheres with slightly different rotation speeds, 
boundary layers of size . are expected near the inner sphere and the 

cylinder circumscribing it (see figure^. 

Note that equations (|1.3j) are linear. Generally, the derivation of the boundary layers does 
not involve the nonlinearity of it only matters in stability questions. In short, one 

can say that the nonlinear term does not create the boundary layer, but may destabilize it. 
Therefore, we restrict ourselves to linear equations 

A'{x,D,)u' = f. (1.4) 

Classically, the problem of the derivation is tackled through one of the following ways: 

i) an explicit calculation, where the exact solution is computed and expanded. However, such 

technique is restrictive (an analytical computation is rarely tractable) and often tedious 
(see for instance with Bessel functions). 

ii) the so-called "method of matched asymptotics" , see 122| I23j for a general presentation. 

The idea of this method is to patch two asymptotic expansions: an "outer" regular 
one, far from the boundary, and an "inner" singular one, close to it. Such expansions 
must coincide in an overlap domain, which provides boundary conditions for both outer 
and inner profiles. The structure of the inner expansion (including the layer sizes) 
is sought by trial, through the principle of least degeneracy (c./. [SI for a detailed 
explanation). This method has been applied with success to various physical systems, 
including singularities near edges. However, the determination of layer sizes leads to 
heavy computations, and must often be supplied with refined physical arguments (see 
|21j on rotating fluids). 

A new method of derivation has been recently introduced in [Sj (see also the proceedings l9l ) . 
It is based on the asymptotic analysis of the symbol A^{x,(^) as e goes to zero, for x £ d^l. 
Briefly, up to use local coordinates near x G dO,, one can consider that equation ()1.4|) holds 
in a neighborhood of x = {x', 0) in . Then, the leading idea of the derivation is to carry a 
Fourier-Laplace analysis: boundary layer sizes near x are deduced from modal solutions 

u%x) = exp (if • x) y^ ^ 0, (1.5) 

of the equation with frozen coefficients 

A''{x,di)u'' = 0. 

In other words, one must consider the characteristic manifold of A^, 

a{a^) = {{x,e),a'ix,e)=0} 

where a^{x,£,) = det^^(x,^). Broadly speaking, if ^ = (CiCn) is the dual variable of x = 
{x',Xn) G , boundary layers correspond to (x,f ) in a (a^) that satisfy 

^''^in ^ +00, e ^ 0. 

The size of the boundary layers is then given by \Xm ^'^\~'^ . In most cases, boundary layer 
equations follow, by appropriate rescaling of the symbol. This method has been used in [S] 
at a fully formal level. It has been applied with efficiency to various geophysical systems. 
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including rotating fluids or MHD. It has allowed to recover the main boundary layers of the 
physical literature, with very few and simple algebraic computations (see for all details). 

The aim of the present paper is to provide this formal method with some mathematical 
basis. We will limit ourselves to scalar equations. In this reduced (but still large!) setting, we 
will show how the microlocal analysis of (|1.4jl is linked to boundary layer expansions (|1.2j) . 
We will present a set of conditions (assumptions (HI) to (H5)) that ensures the existence of 
such expansions, with clearly identified boundary layer sizes and equations. Hypothesis (HI) 
to (H5) involve symbols derived from a^. They will be shown to be intrinsic, although these 
symbols depend on the choice of local coordinates. 

Let us specify the framework of the study. Let X be a smooth Riemannian manifold, 
n = dimX. In all situations of physical interest, X will be of the type T"i x M^^, ni + n2 = n, 
where the torus T"i models periodic boundary conditions in ni physical variables. The 
Riemannian metric will be the Euclidean metric induced by ]R^"i x R"'^. Let Q, a domain of 
X. We denote Y = 90 the boundary of O, which is a n — 1 dimensional submanifold of X. 
We assume that Y has a finite number of connected components, Yi to Yat. We endow Y with 
the Riemannian structure induced by X. We assume that there exist a tubular neighborhood 
of y in X, T > 0, and a smooth diffeomorphism 

T^Yx{-T,T), x^{y,t), 

such that t > in Tnri, t < in Tn (X\Jl). For instance, such assumption is satisfied when 
X = and 17 is a bounded open subset. In the whole sequel, we identify T and Y x (—T,T). 
We consider the scalar problem, 

a''{x,D^)u^ = f, x£n, (L6) 
bl{x,D,)u' = gi, xeY, /G {!,..., L}, (L7) 

where a^{x,Dx), respectively bi{x,Dx), is a differential operator depending smoothly on 
X G n U T, respectively x £ Y, and polynomial in e. The source term /, respectively gi, 
is smooth on QuT, resp. on Y. Let m = deg(a^) > 0. 

In the region Y x [0, T'), T' > small enough, we can write 

m 

a'{x,D^) = a%x,Dy,Dt) = ^a^{x,Dy)Dl 

j=0 

where is a differential operator on Y, of order m — j, with coefficients smooth in x, 
polynomial in e. We assume that the leading coefficient reads 

a^(x) = e*^a^(x), M G N, a„(x) / 0, Vx. 

Before we state our main result, we still need to precise what we mean by boundary layer 
expansion. Let 99 be a cut-off near the boundary. Precisely, ip G C°°{Q), (p = 1 on y x [0, T/4), 
and = on 17\ (y X [0,r/2)). 

Definition 1. A family of functions u^{x) is of boundary layer type if it reads 

u%x) = <(x) + ip{x)vli{x) (L8) 

where the regular part and the singular part have the following asymptotic expansions, 
for some positive 6 and 7* ; 

00 

u'^{x) ^ ^e^^u^{x), u^eC^fpi, (L9) 

fc=0 
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uniformly on every compact subset of $7, and and 

vli{y,t) ~ Y.e'>^ l^'(y)^f (^^'4)' (1-10) 

fc=0 l<i<Af ^ £ ^ 

i<i<n 

uniformly on every compact subset ofYx [0, T), where Vj'^ = v^j^{y,6) are not all zero and 
belong to X' = S [W^; C°°{Y')) . 

We recall that 5 stands for the Schwartz space of fast decreasing functions. Spaces x* 
express the localization of boundary layer profiles. Note that in our definition, the singular 
part is non zero (at least one u*'^ is non zero). Our main result resumes to 

Theorem 1. Suppose that assumptions (HI) to (H5), c.f. next sections, hold. Then, the 
system 

a%x,D^)u'' = f, xen, 
bf{x,D,)u' = gi, xeY, le{l,...,L}, 

has a solution of boundary layer type, modulo e°° uniformly on every compact subset. 

Theorem n is the rigorous translation of the formal method described in article |H]. In the 
following, we will present assumptions (HI) to (H5), which allow to build solutions of bound- 
ary layer type. Broadly, assumptions (HI) to (H3) (see section relate to the determination 
of exponents 7j. Assumption (H4) (section is linked to the derivation of boundary layer 

profiles v^j^. Finally, (H5) (section^ is connected to the construction of the regular part 
u^. The proof of theorem ^follows. Application to the quasigeostrophic equation is given in 
section [3 



2 The singular perturbation 

2.1 Local characteristic manifolds and assumptions (H1)-(H2) 

The aim of the paper is to describe the structure of the solutions of near the boundary. 
Let y £ Y, and (O, x = i^' ^ ^n)) a local chart in X around y, with 

xniy) = o, ^(y)/o. 

The local coordinates {x',Xn) lie in M"^^ x M. We denote by = a^(x,^) the symbol of the 
operator in this local chart: precisely, 

a^^(x,0 := e-*>^(")-« 

Our idea is to deduce the singular structure of from the local symbols a^, precisely from 
their characteristic manifolds. We are interested in modal solutions that are singular with 
respect to e in the direction normal to the boundary. This means we wish to consider wavevec- 
tors = iCi^n) such that 

Id ^ +00, e ^ 0. 

We have the following 
Proposition 1. For all C £ M"~-^, 



a%x,Di)e'^^^>^ 
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i) the m roots of a^(y, C) ') can be written as m functions (e), • • Cm(e) for e > small 

enough. 

ii) There exists p = p{C,) G N* such that, for i = 1, . . . ,m, (e^) has an extension meromor- 

phic in £. 

This proposition follows from standard results of complex analysis. For fixed a^{y,^,') 
is a polynomial in ^„ with coefficients holomorphic in e, and we refer to Kato 13^ for a detailed 
study. 

Corollary 1. For all ( G R""-*^, 1 < i < m, one of the two following possibilities occurs: 

1. ^iie) has a limit as £ > goes to zero, that we denote ^^q. 

2. there exists a unique pi = Pi{0 G Q^!' such that e^* ^i'(e) has a non-zero limit as £ > 
goes to zero, that we denote r/^Q . 

Note that if the leading coefficient of a^(y, •) does not vanish when e — > 0, the second 
possibility does not occur. Roots with singular behaviour exist only if the highest order term 
goes to zero with e, which is typical of a singular perturbation. 

We denote by s = the number of roots satisfying 2. Several indices i may cor- 
respond to the same value of pi. Let 71, . . . , 7,., r < s, be the distinct values of pi. For 
1 < J < we call rrij the number of indexes i such that pi = 'jj. We make the following 
assumption: 

(HI) There exists R = R{y) > 0, such that for all \C, \ > R, the values of r, (71, mi) 
up to {jn^r) are independent of 

This also implies that s = '^rrij is independent on Q. Note that, using notation 7^,, we 
implicitly assume r > 1. When r = 0, the statements in the sequel become empty, so that we 
do not pay attention to this case. 

Assumption (HI) allows to state 

Proposition 2. Let \C,\ > R. There exists W = W (C) « neighborhood of C, O'^d constants 
^0 = eo (C) > 0; = (C) > 0, such that: for all C in W, for all < e < Eq, the roots of 
a^{y, C, •) can be divided into r + 1 disjoint sets Zj{Q, < j < r satisfying: 

i) card 2^o(C) = rn — s, and for all ^„ in Zq{C,), 

< ICnl < 5-\ 

ii) For all \ < j < r , card 2^jiC) = i^ji CLnd for all in Zj{(^), 

Remark 1. Suppose proposition\^ holds. Let (" G W. By propositions ^\ andQJ for e small 
enough, the roots of a^^{y,C,i') can be written Cj^(e), 1 < i < m. Up to reindex the roots, they 
satisfy: 

7^ ^Jo' s+l<i<m, 
£.j vio / 0, 1 < j < r, pi = 7j. 
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Therefore, there exists Eq such that , for all < e < e^, 

m) = {ehe), s + l<i<m}, 

for 1 < j < r. We emphasize that in proposition\^ the sets 2^j{C) o-re defined for < e < Eq, 
where Eq depends only on C- This is not the case for the Cj^(e), defined only up to Eq, which 
depends a priori on Moreover, up to consider a smaller eo, sets -2^|(C) satisfying i) and ii) 
will he automatically disjoint. 

Proof: Let \Q\ > R. With the notations of remark ^ for all < e < Sq, we can factorize 

r 

4(y,C,en) = a^o'^(e„) n i^^'^n) , (2.12) 



where 



al'^i^n) = e^'e-^^^^^C^ n - ^H^)) > 



i=s+l 

afiv) = n - ^^'^^(^)) ' for 1 < j < r, 

P»=7j 

~Q^iy)] ^m{y)' We then define 

with f3r = M — m^r- It is a polynomial in ry, whose coefficients are functions of the type 

ko 



J2 ^'^(o 

k=—ko 



ak 
E , 



for some fixed /cq ^ a > 0, and polynomials symbols 5^. 

Using (|2.12|1 . it is easily seen that: for all locally uniformly in 

It is a priori pointwise convergence in ^. However, the pointwise convergence of the coeffi- 
cients, of the type 

J2 UC)e"' — > 6(C) 

fc=— fco 

holds if and only if 6^ = for A; < 0, and 6o = b. This implies that the coefficients of 

Pi=lr 

are polynomials, and that the convergence is in fact locally uniform in C,. 
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Let \C\ > R. Let Tr a curve enclosing the rj^Q and no other root of a^.^^{C,-). As this 
polynomial has smooth coefficients, Rouche's theorem yields a neighborhood W of C such 
that, for all ( in W, Tr encloses the ry^Q (and no other root of a^^^iC, •) )• Now, the convergence 
of ar,x ^r,x locally uniform in (. Still by Rouche's theorem, and up to take smaller W, 
we get: for all C in W, encloses exactly roots of ar,x(C) ')• 

Back to a^, this yields a 5 = 5{C) > 0, such that a^(y, •) has rrir roots satisfying 

This provides us with the set Z^{Q. 

The construction of the other sets is made inductively, using the polynomials 



with Pj = M — m'jj + Ylk>j+ii'yj ~ Ik)- We only indicate how to build Z^_-^{C). The general 
induction argument is left to the reader. 

The coefficients of the polynomial a^_;^^(Cr) a^^e of the same type as a^^x^Ci')- Using 
again (|2.12)) . we have this time, locally uniformly in Cr/, 

iCv) - 

a'r.,jc,v) := n - io) ( n -io) , 

Pi=7r— 1 \Pi=7r / 

and a^-i x smooth symbols in as coefhcients. Note that 

n -io + 

Pi=7r 

by definition of the t/^q. Reasoning as above with Rouche's theorem, up to reduce W, Eq and 
^, we find a curve F^-i (independent on in W, e < eo) which encloses exactly rrir-x roots 
of a^(y, C, •), all satisfying 

< \in\ < — -, C G W. 



e7r-l — — e7r-l 

They define the set Z^_i{C)- □ 

Let Ci VV and Eq as in proposition |2 For all C in VV, for all < e < EQ) we can factorize 
the symbol as 



(y, C, en) = a§(C, Cn) n «i ('^' ^^'■^") ' (2.13) 



where 



a'oiCU) := C'x n - 

«5(c>^) := n - ^^'■^)' 
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for I < j < r. Note that for all C in VV, and for < e < Eq, we have 



'0' 

where a^''' was introduced in (|2.12|) . We can extend to e = by 

m 

i=s+l Pi=lj 

for 1 < J < r. We emphasize that all are globally defined in Q, i.e. on {|(^| > i?}. On the 
contrary, the are only locally defined in (e, C); on an open subset [0, eo(C)) ^(0- The 
regularity of these various polynomials is given in 

Proposition 3. For all < j < r: 

i) The coefficients of are symbols in (restricted to \(\ > R). 

ii) Let \C,\ > R. The coefficients of a^j are smooth functions of (e, C) o, neighborhood of 

{e = 0,O. 

Proof: We use the notations introduced in the proof of proposition |2j Again, we deal 
only with j = r,r — 1, and leave the general induction argument to the reader. 

i) In the course of the previous proof, we have shown that the coefficients of 

are polynomial in ^, |^| > R, so that the coefficients of aj! share the same property. Also, the 
coefficients of 

<'?-l,,(C.'!) = C, ,•»-""-•»' 



■ n (" - "Jo) i n -4) 

Pi=7r-1 \Pi=7r / 

C^r^"-""-^ a'r-iiCv) ( n -ioj ■ 

\pi=lr / 



are polynomial in (. As 

a 



X 



n -"f,o) 

0i=7r / 



is the coefficient of order m — nrir of ar,x' ^ smooth symbol, which does not cancel by 
definition of the t/^^q. We deduce that the coefficients of a^_^(C, •) are smooth symbols in C, as 
well. 

ii) In the course of previous proof, we have shown the existence of a curve T^- (independent 
on C G W, e < So)) enclosing the set e^'' Z^{C)-, and none of the other roots of ar,x(^' ")■ 
Cauchy's formula, we deduce for all A: G N, 

The right-hand side, so the left-hand side, defines a smooth function of (e, C) in [0, eo) ^ 
Now, it is well known that the coefficients of a polynomial are themselves polynomial in 
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such symmetric functions of the roots. We deduce that the coefficients of a^(C, •) are equahy 
smooth. We can proceed similarly for a^_^, through the formula 

and the smoothness of o^-i follows as well. □ 

On the basis of (HI), we hope for an asymptotics of type (|1.2j) with aj = e^^ for some 
indices j. Nevertheless, for such an asymptotics to be true, it is reasonable that r, (71, mi), 
. . . , (7r.,mr) depend neither on the local chart X) nor on y G y. With above notations, we 
have : for all |C| > -R, for all 1 < j < r, 

a°(C,0) = (-1)-^ n + (2-14) 

By proposition 131 a^(C,0) is a symbol in (" (restricted to \C\ > R). By equation I2.14[ it 
does not cancel. In general, this property of no cancellation is not preserved by a change of 
variable. Therefore, the assumption (HI), as well as r,(7i,mi), . . . , (7r,mr) depend on the 
local charts. To overcome this problem, we need to make a stronger assumption: 

(H2) For all 1 < j < r, a°(-,0) is elliptic. 

Note that, up to consider a larger R, (H2) implies (|2.14j) . Note also that (H2) is preserved 
by a change of the tangential variable x'. 

2.2 Asymptotic invariance and assumption (H3) 

We can now state the following invariance result: 

Theorem 2. (Invariance through diffeomorphism) 

Assume (H1)-(H2). Let {Ox^,'^) a local chart in X around y such that 

M/„(y) = o, ^(y)/o. 

Let (H1)^-(H2)^ he the same as (H1)-(H2), with in place of x- Then, (H1)^-(H2)^ holds, 
with the same r, (71, mi), (7r,m,.). 

Remark 2. In short, theorem\^ states that the conjunction (H1)-(H2) is intrinsic, as well as 
the main features of the singularities (exponents 7^ and multiplicities r,mj). We stress that 
(H1)-(H2) involves the whole symbol a^{y, •) and not only its principal symbol. Therefore, it 
is not obvious that it should be preserved by another choice of coordinates. 

Proof: Let us first set a notation: for any smooth (p = ((/>', (pn) : U ^ V, €z U, for any 
symbol P{C, rj) defined on M'^-i X M, and for any 7 G M, we introduce 

We also denote 

Let {0,x) and (0\j,,^') the local charts of the theorem. It is clear that (HI) and (H2) 
(and r, (71, mi), . . . , (7r, m^)) are not affected by a translation of each coordinate. Therefore, 
with no loss of generality, we can assume that xiu) = ^(y) = 0- Let = o x~^- It is a 
smooth diffeomorphism between two neighborhoods U, V of 0. Let now P^{C, v) a polynomial 
in e, C, rj. The proof of theorem |21 will follow from the study of (fiP^- Precisely, we will show: 



1 

2iTT 
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Proposition 4. Let 7 > 0. Then, uniformly in each compact subset of . 



pn— 1 



6^0 



(2.15) 



where (f){x' , Xn) = ((/>'(x', 0), 0) x„) . 

Remark 3. If we denote 

the expansion of with respect to t], then 



with 



(2.16) 



\x'=0 



The symbols are polynomial in More precisely, classical symbolic calculus (see 0) 
provides the following (finite) expansion: 



a>0 



Jr'ix') [!^(^'^Qy 



dXn 



\x'=0 



D^P^ (d.'</.'(0)C) 



(2.17) 



where 



r' {x') := (l)'{x',0) - (i^/0'(O) {x') . 



The terms of rank a in this expansion are polynomial in (, of degree less than deg yPjj ~ 10^1/2- 
(see again for details). 

Remark 4. In the case 7 = 0, the limit ()2.15j) is easily replaced by 

lim,/.,P^(C,r?) = cP.P'^iCr]). (2.18) 
e^O 

We postpone temporarily the proof of the proposition, and we show how to deduce theorem 
121 from it. We must prove that (Hl)ij,-(H2)^ holds, that is we can replace the local symbol 

by the other local symbol a|,. With a slight abuse of notation, we identify x and x{x) so 
that we work with x £ U, U a neighborhood of y = in M". Hence, 



(2.19) 



Proof 0/ (Hl),j,: We consider 



and, as in previous proofs, 
Equation 1)2. 19(1 becomes 



a^_^(C,r/) :=e-^'-a|(0,C,^), 



<,^(C,r?) :=£-^^a^(0,C,^). 



We have seen in the previous proofs that a^^^ is polynomial in e",C,r/ for some appropriate 
a > 0. Moreover, 

aO^(C,7?) = C^??"aO(C,77), |CI > i?- (2-20) 
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If we expand 

j 

we get from proposition 0] that, for |C| large enough, 

hm 4,^{C,v) = ^b'^rjiO rf , 



with 



x'=0 



Using (|2.2()|) . we deduce that the term of lowest degree (with respect to t]) in lim a^,j,(C,ry) 

e— >0 ' 

is b^m-mriO ry™'"™''. Moreover, by (|2.17j) . the principal symbol of b^m-nir 

= C^a{a%,0)) (#^,(0)C). 

Note that dx'4>'{0) is a diffeomorphism in R"~^. By assumption (H2), aj.{-, 0) is elliptic, hence 
brm-mr ^^^^^ elliptic. As a consequence, it does not vanish for |^| large enough. Therefore, at 
fixed C, the polynomial lime^o o,r^iCy ') has rrir roots away from zero. By Rouche's theorem, 
it is still the case for e < Eq small enough. Back to the original variables, we find rrir roots of 
a|(0,C,-) satisfying 

4 ^ < '4- 

With similar arguments, we can find nir^i roots of a|,(0, C, •) which scale like e""'^^'^. We 
consider this time 

"r-i,* •= ^ a|,(0, C, )' 
as well as a^_i -^. Again, a^-i x ^® ^ polynomial in (e",C,ry) and 

a^r-iJCv) = C^r?—— '-aO_i(C,r?)aO(C,r/). 

We then apply the same reasoning as above, replacing a^(-, 0) by the product a^_i{-, 0)o[l(-, 0), 
which is still elliptic. 

Proceeding recursively, we find, for all 1 < j < r, nij roots of a|,(0, C, •) such that 

^ < u < 

for e < eo small enough, 6 = S{eo) small enough. Finally, by (|2.18l) . 



hni a%{0,C,U) = <t>* I aliCU) JJa^jiCO) 



At fixed the right-hand side defines a polynomial in of degree m — s. By Rouche's 
theorem, this yields m — s roots of a|,(0, C, •) with 

< \U < 6-\ 

As the degree of a|,(0, C,-) is m = m — s + Yl''^ji we obtain in this way all the roots. 
Assumption (Hl),i, follows with the same r, 71,. . . , 7r. 
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Proof o/(H2),i,: We have just established (Hl),i,. Just as for a^, we can factorize 

r 

aUyX,U) = aoCC^n) Ha'jiCe'y^U) (2.21) 

i=i 

(see factorization ()2.13|) ). and we have to show that: for all 1 < j < r, a^{(^,0) is elliptic in a 
neighborhood of y. But from above identities, we have got easily: 

C.i,a(aO(-,0))(C) = C^a(aO(.,0)) (#^,(0)C), 

and 

^7(aO_i(-,0))(C) ^T(aO(.,0)) (C) = 

a(aO_i(-,0)) (4'</.'(0)C) a(a°(-,0)) (d.'</)'(0) C), 

and so on. As ellipticity is preserved by inversion and product, we deduce (H2),i, from a 
simple recursion. □ 

Proof of proposition [l] The proof relies on the representation of (pJP^ as an oscillatory 
integral. The asymptotic analysis of such integral will be performed through "stationary 
phase type" theorems. This approach is very classical in symbolic calculus, to establish the 
properties of conjugation, product, or transformation under a diffeomorphism (see, among 
many, textbooks [SlOl)- However, in the standard setting, one is mainly concerned with the 
principal symbol, so that the small parameter in the asymptotics is the wavelength 1/|^|. In 
our framework, the natural idea is to work with the parameter e instead of the wavelength. 

Let ueC^iU), C ^ r/ / 0. We wish to compute 



\x=0 



P = [p^(D,.,e^Z?,j(e-^'^("HC,.-^)„(^)y 
P = P{C,,rj) is defined by the oscillatory integral: 

Following the scheme of [H page ??], we make the change of variables: 

dx„(l)nix',0) 

rj ■= rj — V- 

By the Fubini theorem for oscillatory integrals, we can write 

1 



(27r)" 

with 



P^ (^C,£'^V + dx„(pnix',0)v^ dxndfj, 

where 

r^(x) := T] (<?!>„ (x) - dx„(j)nix' , 0)) . 
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We set 

^%x',fj) := J e"""*' u{x) e^"^'(^)< dx„, 

so that 

r{x',0 = j (C, e^?? + 9,>„(x',0)r?) $^(x',ry)dr?. 
For fixed r/, x', we introduce the phase 

^''^ f{xn) ■■= e~'^ (rnix) - e^'7)xn), 

with derivative 

f{Xn) = V {dxn<Pn{x) - d^^4>n{x', 0)) - E^'fj. 



• li\r] dx„(pn{x' ,0) + e'^fjl < 1/C, C such that | r/(?a;^(?!)(x) | > 2/C for x in a neighborhood 
of 0, then \ f'{xn)\ > l/C. 

• If I r] dx^cpnix' , 0) + e^fj \ > C, C such that | r]dx„(j){x) \ < C/2 for x in a neighborhood 
of 0, then I f'{xn) I > C/2. 

We deduce from the non-stationary phase theorem (see [21 page ??]) that for large enough 
C, outside 

^ < \vdxMx',0)+e^v\ < C, 

for all k, 

\^%x',fj)\ < Cfe + + (2.22) 

Still following [2j, we introduce a truncation function a such that a S C^(M), a(.^„) = 1 
for C7-1 < |e„| < C, a(e„) = for \U\ < {2C)-^. We divide 

.rix',c) = .mx',o + Jiix',o 

with 

J| = ja (^dxMx',0) + ^r/^ (x, C, e^r/ + 9,,>„(x', x„) r/) $"(x',r?) dr?. 

Note that and are polynomials in C (as the only dependence on ( is in P^). From the 
estimate (|2.22|) . we deduce 



Hence, 



If := j^le-'<J!{x'X)dx'dC 

satisfies If = 0{e^) for all k. Moreover, a rapid look at the proof shows that all estimates 
hold uniformly in every compact subset of M"""*^ x M*. 

It remains to estimate the "stationary terms" Jf and If. For this, we will use a lemma 
whose statement and proof can be found in |2j: 
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Lemma 1. (l^) 

Let r = r{z,d) e C°°{R x O), O open subset of W , such that d^r{0,9) = 0. Let = 
X M X O), depending on a parameter X > so that, for some M , 

\d^d^'b\z,v,e,X)\ < Cp^p, X''-"'. (2.23) 

Assume moreover that b^ has compact support in z, uniformly with respect to rj, 6, X. Then, 
the integral 

J(A) = (27r)"" f e-'''^e'^''^''^H^{z,r,,e)dzdr] 



has the following asymptotic expansion for X — > +oo 



\z=0, 



a>0 



where terms of rank a are bounded by A*^""/2. This expansion holds locally uniformly with 
respect to 9. 

The function Jfl^iC) is polynomial in C,, with coefficients that read 

-ix„ f, gie-^r,(x) (-^ 

Such integrals can be easily put in the framework of lemmaQ with A := e^'^ , z := Xn, rj := fj, 
and := {(,r],x'). Thus, r{z,9) := rr,{x) satisfies dzr{0,9) = 0. Besides, the function 

b^{z, 7], 9) := WAxn, ??, C, x) 



satisfies estimate (|2.23|) with M = 0. Its support in z is contained in the support of u, so 
uniform with respect to the other variables. Hence, we can apply lemma ^ The first term of 
the expansion provides 

J'2{xX)^4{xX) ■■= e*^'(-'''°)P°(c,5,>n(x',0)r?) u{x',0) 

and consequently 

II ^ j^le-'<4{x,C)dx'dC 

with convergence in every compact subset of R"^^ x M.^. As we handle polynomials in rj, the 
uniform convergence extends to compact subsets of M"^^ x M. By taking u equal to 1 in a 
neighborhood of in U, we obtain the result. □ 

We end this section on asymptotic invariance with a definition: 

Definition 2. Let y €zY. We say that a^ is uniformly singular at y if (H1)-(H2) holds at y. 
The corresponding jj 's, 1 < j < r are called singular exponents (of a^ ) at y, and the ruj 's 
are their multiplicity. 

These definitions make sense because of theorem |21 We recall that the singular exponents 
are positive rational numbers (see corollary^. 

Definition 3. Let Y' C Y . We say that is uniformly singular on Y' , if it is uniformly 
singular at all y C Y' , with constant singular exponents and multiplicities. 

In the whole sequel, we will make the following assumption: 

(H3) The operator a"^ is uniformly singular on each connected component of the 
boundary Y. 
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3 Boundary layer sizes and equations 
3.1 Definitions and hypothesis (H4) 

Let Y' a connected component of Y. By assumption (H3), is uniformly singular on Y'. 
We denote by 71 to 7,. the singular exponents, and mi to 771^ their respective multiplicity. In 
a region Y' X (0,r'), equation reads 

a''{x,Dy,Dt)u^ = 0. 

Let y G Y'. Let {0',x') a local chart in Y' around y. The application 

x:0'x(-r',r')^M", iy',t)^ix'iy'),t) 

defines a local chart in X around y. Hence, we can define a local operator a^{x, D^' , Dt), 
with symbol a^{x, (, ^„). 

On one hand, we have a finite expansion of the type 

a^LDy,^)= X; Y.^>',k'iy,Dy)D';e'^' (3.24) 

^ ^ k=-ko k' X' 

for some appropriate a > and smooth differential operators on Y'. Hence, 

a^y,C,^) = E (3-25) 
fc=— fco fc' 

On the other hand, with notations introduced in l|2.13|) . we have for |C| > R{y), 

lim e^M^ (y, C, ^) = r?— S.'>. ^(C, v) U 4(C, 0) (3.26) 

Prom (|3.25j) and (|3.2(i|) we deduce 

lime/^^a^(y,D„|^) = d';^^^'^^"''^ a,^{y, Dy, De), (3.27) 
where a^^ is a smooth differential operator. It satisfies, for |C| > R{y), 

j'>j 

For 1 < J < r, we say that the operator a^. is a singular operator (associated to the singular 
exponent 7j). 

Special attention will be paid to localized solutions of equation 

a^u = f, (3.28) 



u*''^ in expansion (|l.lfl|) will satisfy equations of type (|3.28|) . Therefore, we introduce the 
Frechet space 



where is a singular operator, associated to the singular exponent 7. Indeed, the functions 

fy equations of type (jSZ 

X := S (M+; C°°(y')) • 

We also denote 

Xq := {u G X, a^u = 0} . 

We set the following definition 
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Definition 4. We say that j is a boundary layer exponent (and a boundary layer size) on 
Y' , if there exist subspaces Xq 7^ {0}, X' with: 

i) Xl^ C Xq, X'^dX' ^ X. 

ii) DeX' = A", DyX' C X' , 0X' C X', f X' C X' for all smooth functions f G C°^{Y'). 

iii) For all f £ X' , (|3.28|) has a solution in X' . 

Remark 5. Note that identity DqX' = X' implies that X' is preserved by both Dg and its 
inverse Dg^ (which is obviously defined on X). Note also that Xq and X' are not uniquely 
determined. However, the subspaces + X^ and + X' where the algebraic sum is taken over all 
subspaces Xq and X' satisfying i)-iii), are well-defined, and are the largest subspaces satisfying 
i)-iii) 

Remark 6. The definition of a boundary layer exponent is at first sight dependent on the 
diffeomorphism 

X-.r ^Y' x{-T',T'), x^{y,t), 
where T' is a tubular neighborhood ofY'. However, if 

i;:f'^Yx{-f',f'), x^{y,i) 

is another such diffeomorphism, we can link the singular operators a^^^ and a^^^ following 
computations of theorem\^ Indeed, we have 

aj,4,iy,Dy,Dt) = '4).^a^^y,{y,Dy,Di), 

where 4> = (</>', (/>„)= -0 o X~'^ £Y' x {-Tmin, Tmin), where Tmin ■= min(r', T"). , and 

Thus, solutions of (|3.28)) with a-y = a^^x ^'^'^ '^7 ~ ^i^i> ^'''^ deduced from each other by the 
change of variable (f) . As soon as 

this change of variables preserves the properties of a boundary layer exponent. 

We end this section with an assumption that will ensure the existence of a singular (bound- 
ary layer part) in p.8jl : 

(H4) Among the singular exponents given by (H3), there exists at least one 
boundary layer exponent. 

3.2 Solvability 

The verification of (H4) relies on solutions of (|3.28|) in X. In some particular cases, one can 
give some effective criteria that ensure the solvability of this equation. 
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3.2.1 Order zero operators 

For many physical systems, the coefficients of the singular operators are just smooth functions 

a-y{y,Dy,D0) = a^{y,De). 

Hence, (|(-i.28|) resumes to a collection of linear ODE's with constant coefficients, indexed by y 
in Y' . Thus, to determine if 7 is a boundary layer exponent is relatively easy. For instance. 

Proposition 5. Assume that there exists a closed curve 

V C {Im [t]) > 0} 

enclosing the set 

{r/, 3y, a^{y, r/) = 0} n {Im {r]) > 0} 

of all roots with positive imaginary part. If this set is not empty, then ^ is a boundary layer 
exponent, otherwise it is not. 

Proof: By Rouche's theorem and smoothness of a^, the number m'^[y) of roots of a-y{y, •) 
lying inside F is locally constant, As Y' is connected, it is constant. It is then well-known 
(see P) that the homogeneous space Xq is spanned by functions of the type 

f{y)i.Oj{y,e), feC^iY'),u;,iy,e) = ^ e^"^ a-i(y, r?) ^ dr?, j G N. 

Moreover, the dimension of Vect {u)j{y,9), j £ N} is m+. The result follows, using Xq = Xq, 
and 

X' = Vect f^fiy) j^e'"^'a-\y,r])rfdr], / G C°°(y'), ^ G N,i E z| 

□ 

As shown in proposition El the boundary layer is connected to roots with positive imag- 
inary part. This echoes the formal method used in where roots .^^ with Im +00 
were considered. 



3.2.2 Operators with constant coefficients 

Suppose that Y' = M"~^ (or T"^^) and has constant coefficients. One can use the Fourier 
transform over M"~^ (or T"~^), which turns ()3.28() into 

a-y{C,Dg)u = f. 

This is still a collection of linear ODE's, with constant coefficients, indexed by ^ G M"^^ (or 



3.2.3 General case and link with spectral problems 

In general the singular operator a-y{y, Dy, Dg) will be a "true" differential operator. But in 
many cases it will give rise to a (spectral) decomposition which will allow us to solve the 
problem. In this paragraph we will describe briefly the methodology without focusing on 
technical details. 

Let us remind that a solution of a problem of type 

a^{y,Dy,De)uiy,e) = f{y,d) (3.29) 
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that decreases as — > +00 can be handled by Fourier-Laplace transform in 9. Denoting 
u{y,^) := e~^^u{y,9)d9, H3.29() is equivalent to an equation of the form : 

a^{y, Dy, u{y, i) = f{y, + F{y, := G{y, (3.30) 

where F involves boundary terms issued from integrations by parts in the Laplace transform. 
Let us suppose that := a^{y, Dy,S^) admits a ( sometimes spectral) decomposition weakly 
on of the form: 

[ g{X,0\^i,x X ^i,x\dX (3.31) 



where \^(,x >< ^(,x\ denotes the orthogonal projector on f^^Xi and the (spectral) parameter 
A is either complex or real. Let us suppose moreover that the function g{\,S,) and the family 
of vectors ip^^x are analytic in ^ (let us mention here that a naive definition of analytic vectors 
is enough for the examples we have in mind: for a precise definition see |17|). We can solve 
(|3.3()|) by using the decomposition (|3.31|) . Indeed denoting 

G(A,0 ■.= <n,x\G{X,.) >:= j ^G(A,e)de 

we get: 

uiy,0) = — d\ d^——G{X,Oe^\^,x{y)- (3.32) 

Therefore by the Cauchy theorem it is enough to know the zeros of g. Let us show how these 
ideas apply in the case of a non-differential operator, a pure differential one and in general 
when the boundary is a flat manifold.. 

Let a^{y,Dy,Dg) = b{y Dg). In this case if^^xiv) = ^{y - A) and g{X,^) = b{X,C). 
Therefore the l.h.s of (|3.32j) can be computed easily: let Ciiv) a family of non-degenerate 
zeros of b{y,C), that is b{y,^i{y)) = 0, d^b{yCi{y)) / 0. then: 



/ dC-—-GiX,0e^'6{y-X) 



!• j-X+lOO Y 

-I f'X-\-ioc 1 

= 17- / dCTr^Ay,Oe'^' 

^■^i Jx-ioo o(y,0 
= E^TT^ G{yMy)y^'^'^' 

One sees immediatly on this example that the boundary layer will exist as soon as one of the 
E,i{y) satisfies TZe{^i{y)) > c > 0, and there is no crossing between the different branches. 

Let us suppose now that a^{y, Dy, Dg) = c{Dy,Dg). In this case ip^^xiy) = e^^^ and 
g{X,^) = c(A,^). Let as before ^i(A) a family of non-degenerate zeros of 6(A,^), that is 
6(A,e^(A)) = 0, a^6(A,Ci(A)) ^ 0. We get: 



dX / dC——G{X,C)e^'e'^y 



27ri 

EtHTT^ G(A,6(A))e 



^^(X)e^iXy 
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Therefore the boundary layer will exist as soon as there exists a family ^i(A) with TZe{S,i{X)) > 
c> 0, VA. 

In the case where the boundary is M" (or T") and that the operator a^{y, Dy, Dg) is a 
differential operator in y with polynomial coefficients a decomposition similar to the preceding 
one can be handled with the so-called coherent states. Let us briefly describe the construction 
on M, the generalization to M" (and T" by periodization) is straightforward. Let (p, q) £ M? 
and let (pp^q the family of functions defined in ^ 

Let us now consider the operator a^{y, Dy,£) = a{y, Dy, ^). Being with polynomial coefficients 
it can be surely rewritten as 

a{y, Dy,0 = + dy)Hy " dyY (3.33) 

Let us define aAw{q,P,0 '■= Efcj + w)''{y - 

Lemma 2. 

aiy,Dy,0 = / aAwiq,P,(.)\Vp,q >< '^p,q\dpdq 
where the integral is understood in the weak sense. 

Proof: an easy computation shows that we have the following decomposition of identity 

Y \'Pp,q >< ^p,q\dpdq = 1 

Moreover it is easy to check that 

{y + dy)Lpp^q = {q + ip)^p,q 

and therefore 

< 'Pp,q\{y - dy) = {q- ip) < iPp,q\- 

Inserting in (|3.33|) the decomposition of the identity we get: 

{y,Dy,o = Y.'"'AO{y + dy)Hy-dyy 

= Y] + 9y)^ / \^p,q >< Vp,q\dpdq{y - dyY 

= / y2o!kj{q + ip)''\(Pp,q><iPp.q\{q-ipydpdq 

aAw{q,P,C)\^p,q >< ^p,q\dpdq. 

It is easy to check that the same argument as for the preceding case holds, and that the 
boundary layer exists as soon as it exists a family of non-degenerate zeros, OAwiq^P, Ciq^P)) = 
with real part strictly positive (uniformly in {q,p))- A detailed description of this method 
will be given elsewhere. 
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3.2.4 Almost solvability 



Contrary to what previous cases suggest, the solvability of (|3.28|) can not be deduced in 
general from a symbolic analysis. Indeed, to go from a local analysis (with a local symbol 
near y E Y') to a global one (on all Y') requires a "patching process", involving a covering 
of the boundary and a partition of unity. Such process introduces regularizing operators, 
responsible for a loss of information at low frequencies. One can clarify this idea in the light 
of classical results by Treves jSUl Chapitre 3]. These results were obtained to describe the 
regularity properties of elliptic operators. However, they give some insight into our boundary 
layer problem. Let us consider the "tangential principal symbol" 

m 

A^{y,Dy,De) := Y,^M{y,Dy)Dl 

where := XlfcLo O'kiy, Dy) Dg. We suppose that 

i) Y' is compact 

ii) There exists q > 0, such that 

degree cr(afc) = degree a (cm) + q{m — k), Mk. 

ill) There are two integers, rri^ and m~ such that m"*" + m~ = m, m+ > 1, and for all (y, Q) 
in T*Y' — {0}, the polynomial A^[y,C^,r}) with respect to r] has exactly m"*" roots with 
positive imaginary part, m~ roots with negative imaginary part. 

Under these assumptions, we state 

Theorem 3. There exists a family Rk, k = 0, . . . ,m of regularizing operators over Y' such 
that 

m 

dy{y,Dy,Dg) := a^{y,Dy,D0) + '^Rk{y,Dy) Dg 

satisfies: {u S X, d^u = 0} 7^ {0}, and for all f £ X, equation d^u = f has a solution in X. 

Broadly speaking, this theorem says that a singular operator whose "symbol has good 
properties" (assumptions i) and ii)) is "close" to a boundary layer operator, i.e. up to 
additional regularizing operators. However, as these regularizing operators do not vanish, 
one can not conclude to the solvability of (|3.28j) in X. 

Proof: We start noticing that 

m—l 

a^{y,Dy,Dg) = am{y,Dy)D^ + ^ak{y,Dy)D^ 

k=0 

with Qm elliptic. Up to compose by a parametrix of am, which would yield an additional 
R(y, Dy)D'^ term with R regularizing, we can assume that am = 1- 

This simplification, together with ii) and iii) enter the framework studied in [SOI- Treves 
considered ii) with g = 1, but its arguments extend straightforwardly to g G N*. In particular, 
we have the "almost factorization" (see pHI pages 157-162]): 

a^{y,Dy,De) = a~{y,Dy,De) a1^{y,Dy,D0) + R{y,Dy,De), 
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where a± := ^ a±(y, Z)^) I^f, R :=Y,Rk{y, Dy) 

k=0 k=0 

The a^, resp. Rk, are pseudo-differential operators of degree q{m^ — k) resp. — oo. Moreover, 
for all (y, C), the roots of 

4(y,C,r/) := J] (T(a±)(y, C) r/^ 

fc=0 

are the roots of A^{y,C,^-) with ^ imaginary part. Thus, from above decomposition, we 
can rewrite equation (|I'{.2<S|) as a triangular system: 



at/V + Ru = f, 



a^u = V 



(3.34) 



Still following we denote — A^^ the Laplace-Beltrami operator on Y' and set A = (1 
Ay)^/^. Through the change of variables 



Dgv^ = A«t;^+\ j = 1, . . . , m+ - 1, 
Deu^ = A^u-'+\ j = l,...,m- - I, 



(3.35) 



we can write p.34() as a first order system 



DgV - A^v = nu + F, 

DgU - A'U = JV 



(3.36) 



with V = (v^), U = (n^), 

nU = {0,...,A^-"'^Ru^), F= (0,...,Ai-'^''/), JV={0,...,A^-'^~v^). 
The A^ are matricial pseudo-differential operators of order q over Y'. Moreover, they satisfy 

det{r,-a{A^){y,C)) = A^{y,C,v)- 

Hence, we can apply [201 theorem 1.1, page 134] on parametrices for parabolic type oper- 
ators: one can find a function IA^{9), with values in matricial pseudo-differential operators of 
order zero over Y' , such that 

DeU^ ^A^U^ ~ in y X {0 > 0}, U^\e=o = 

where is the identity of C"*. In each local chart (O', x') of Y\ the local symbol (y, Q) ^ 
Ut{y,0,O ofU^ie) satisfies 

i) U^, is C°° on O' x M+ x R"-\ 

ii) To every compact subset /C of C, to every pair of (n— l)-tuples a, /? G Z"^^, and to every 

pair of integers r,N > 0, there is a constant C > such that for all y in /C, > 0, 

\\d^d^^dlUp{y,9X)\\ < C0-^(1 + |C|)'^^-1^'-^^ (3.37) 
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We quote that in theorem 1.1 of Treves, estimate (|m-{7j) holds only locally in 0, because 
the symbol A = A{9) under consideration has only local control in 9. In our case, At^ are 
independent on 9, so that (|3.37j) extends to 9 in M4.. 

Suppose now that f e X so that F G X""^ . Let V° G (C°^(y'))"'^- Set 



1-9 

V = U+{9)V^ - / U+{9 - s)F{s)ds, 
Jo 

r+oo 

U = - U-{s-9)JV{s)ds. 



(3.38) 



By estimates (|3.37|) . it is clear that {U,V) G X"^, and as are parametrices, 

DeV - A+v = nu + n-^v + F, 

DeU - A-U = JV + n-U 
for regularizing operators TZ^. Back to original variables, we get the result. □ 



(3.39) 



4 Boundary layer expansions 
4.1 W.K.B. Ansatz 

The notions and hypothesis introduced above are linked to boundary layer expansions for 
solutions of (|1.6|) - (|1.7() . Let to the connected components of y = dU,. To each 
i = 1, . . . , A^, we can associate G N boundary layer exponents < . . . < 7*,. Remind that 
these exponents are positive rational numbers. Hence, replacing for simplicity e by with 
some appropriate (5 > 0, we can assume that 6 = 1, 7* G N*, Vi,j. 

As usual in WKB approach, we inject approximations H1.8|) - (|1.1U() into (|1.6() . The resulting 
equations are expanded and ordered according to powers of e, and coefficients of the different 
powers of e are set equal to zero. It leads to a family of equations on the and v^-^. To 

lighten notations, we set vf' = 0, v'jf = for A; < 0. 

i) Interior terms. Denoting a^{x,Dx) '■= 'Y^=o^^^i{^i^x)i we have 

h 

ai{x, D,) u^~' = 5ok f, Vfc G N. (4.40) 

1=0 

where 5ok is the Kronecker symbol. This reads in short 

A{x,D,)U = F (4.41) 
where U := {u'') G C°^{Tl f, F := (/,0, . . .), and 

A-.c^inf ^c°^(nf, {Auf ■.= Y,aiu^'^- 

ii) Boundary layer terms. We deduce from 1)3.27(1 

of a^. (y, Dy, De) vf + L(y, 6, Dy, De) (vf, . . . , = (4.42) 

for y £ Y^ and 9 G M+. L is the linear differential operator coming from lower order 
terms, and n*- is an integer. L has smooth coefficients, polynomial in 9. As Dg is 
invertible from X^ to X^, we rewrite such equations as: 

a^.{y,Dy,De)vi^'' + M{y,9,Dy,De) (^^f, • • • , t^f =0. (4.43) 
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In particular, 

a^.iy,Dy,De)vf = (4.44) 

Moreover, as 7*- is a boundary layer exponent, equation H4.44() has a non-trivial solution 
in (^*). Recursively, equation ()4.43|) (hence equation (|4.42|) ) has a solution in (^*)- 
Again, equations H4.42() can be written 

/ N \^ 

Au{y,e,Dy,De)V = 0, V = (y\\ . . . ^vf^) ^ ^ & ' 

Remark 7. No other settlings in boundary layer expansion H1.1U() could lead to "reasonable" 
localized solutions. This means that, in view of factorization (|2.13|) . any boundary layer term 
of the type v{y,t/e^) where 7 is not a singular exponent yields equations of the type 

e{y,Dy)D^v = f 

for some elliptic operator e{y,Dy) and a G N. Except for some spurious source terms in the 
Cttse where e{y, Dy) is non-invertible, such equation has no nontrivial localized solution. 

4.2 Assumption (H5) and proof of Theorem [T] 

At this point, the boundary layer and regular parts have been described independently from 
one another. They are of course linked through the boundary conditions (|1.7|1 . It can be 
written 

m' 

^ Bliiy,Dy)D'su\t=o = 9i, y G Y, (4.45) 

k=0 

where Bf.i{y,Dy) are smooth differential operators, polynomial in e. Let us introduce 

/N sm'+l 

dV = {[DlV\t=o) , k = 0,...,m', AuV = o] C IJlC^iY'fj 
Again, we can inject the Ansatz (|1.8j) - (|l.lfl|l into (|1.7)) . This can be written as 

m 

Y,l3M{y,Dy)D^sU\t=o + Ci{y,Dy)W = d, 

k=0 

for some differential operators Bm : C°°{Yf ^ C^iYf, W e dV, and 

m+l 



Introducing any projector 11/ on the range of Ci, we can rewrite previous equation as 

m 

J2^iBki{y,Dy)D'^U\t=o = liiGi (4.46) 



A:=0 

The existence of the regular part resumes to: 
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(H5) One can find U G C~(J7)^ such that 

(A{x,D,)U = F, 

m 

^UiBki{y,Dy)D^,U\ 



t=Q 



(4.47) 



^ k=0 



As assumptions (HI) to (H5) are fulfilled, one can find of type (|1.8|1 - (|1.1()|) satisfying 



b'i{x,D^)u gi, xGY, l£{l,...,L}, 



(4.48) 



uniformly on the compact subsets. Theorem ^ follows. 

One can see system ()4.47() as the limit system of (|1.6() - ()1.7p . It means that the small-scale 
effects due to boundary layers are contained in this large-scale system. A strong difficulty 
is that it is infinite- dimensional: correction terms it* at any order may interact. Hence, in 
most situations, one does not treat directly system ()4.47|) . The classical attempt is to turn it 
into an infinite collection of finite- dimensional systems, that are of the same type and solved 
recursively. This will be illustrated on the quasigeostrophic equation in the next section. 



5 Application to the quasigeostrophic model 

The stationary linearized quasigeostrophic equation reads, for x in a two-dimensional domain 

n 

Pd^^ip + V-^V-Vr/B + rAip - ^^'^^p = /3curlr (5.49) 

Re 

where ip = ip{x) G M, x = {xi,X2) in cartesian coordinates, f3, r and Re are positive param- 
eters, and T]B = r]B{x), r = r(x) are smooth functions. It is a celebrated basic model for 
oceanic circulation: ^ is a stream function, associated to horizontal velocity field u = SZ-^tp, 
TjB is a bottom topography term. The expression P dxiip comes from the variation of the 
Coriolis force with latitude, rAxjj and — 1/ReA^^ are dissipative terms due resp. to friction 
and viscosity. Finally, /?curl r is a vorticity term created by the wind. We refer to for 
all necessary details. Equation 1)5. 49() is fulfilled with Dirichlet boundary conditions 

^ = = 0, x£ dn. (5.50) 

We will apply previous analysis to system (|5.49|) , (|5.50j) , for several domains Q and ranges of 
parameters /?, r and Re. 

5.1 Munk layers 

We first investigate the case where 

^ +00, r. Re, r given (5.51) 
which corresponds to strong forcing by the wind. We consider a domain of the type 

^ ■= {Xw{x2) <Xi< Xeix2)} , 

where Xw and Xe ai'e smooth, with all derivatives bounded. They describe the western and 
eastern coasts. We assume that X2 G M or X2 G T. The case of a closed basin will be evoked 
in last subsection. 
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This boundary layer problem enters the framework given in section |21 with e = (5 ^ . We 
start with the derivation of singular exponents for both connected components of Y := dVt, 

Yw ■■= {{Xw{x2),X2), X2 G M or T} , y; := {(xe(x2),X2), X2GMorT}. 

Western boundary layers. For any y ^Y^, global coordinates x'^ = X2, x'2 = Xi — Xw{x2) define 
a local chart near y. In these coordinates, the corresponding symbol satisfies as |^i| + IC2I — > 
+00, 

au,{x', ii,i2) := ie-H2 + 0(|6| + |6|) - r (il + (l + X^')^' 

+ + 161)) - ^S^t + i^+x'.?e2 + om\ + Mf) 

with Xw '■— Xwix'i)- For all ^ G M, the roots of a^(x',6, •) with |,^2| +00 satisfy 



1/-. /0\') J.A ■ j~ ( zR,G 



Eastern boundary layers. The derivation is similar, with x'l = X2, X2 = Xeix2) — xi. Roots of 
Oe(x',^i,-) with singular behaviour satisfy 

\{i+x?r) 

Hence, assumptions (HI) to (H3) are fulfilled, with the singular exponent 7 = 1/3. 

To identify singular operators and boundary layer sizes is easy using the global coordinates. 
Indeed, 

Thus, the singular operators are 

a«,,7 ■=^s-^{l + X'^f Dl, ae,^ := -De " (l + x'ef 

They have order zGro coGfficiGnts. Denoting g : — Re^/^/(l + Xwefi the roots of their 
symbols are a^jC*'^/^, a^e^*'^/^, a^jC^*'^/^, resp. aee"""^/^, Oge*'^/^, Oee'^^^^^. Proposition [51 
applies (two roots, resp. one root with positive imaginary part). Thus, e^^^ is a boundary 
layer size on Y^^e- 

It remains to show the existence of solutions similar to ()1.8p . with 

nL,(x) ~/(:e) + e'/^^\x) + ... 

^ ,0 f Xl- Xwix2) \ , ,0 f Xe{x2)-Xl 
Vbl{x) ( ^ ,X2 j + ^Pe ( ^ ,X2 



, .1/3,/.! ( XI - Xw{x2) 

Equation 1)5. 49(1 yields 



X2 + 



d,,i^' = f + 5oiCurlr, (5.52) 



27 



where /* involves , k < i — \. Equation (|5.5()|) yields 



De'4'w\e=o = 9w, ^eV'el0=o = 9l, 

where gf^ g depends on and ^, k < i — 1. The western boundary layer operator has 
two characteristic roots with positive imaginary part. One can find localized solutions ipy^ of 
aw,'yfpw = fw for arbitrary values of 'ipw\e=o^ -C'6»V'«)|e=o- This is not the case for the eastern 
boundary layer, for which only one condition can be given. The last boundary condition 
prescribes DqiP'^Iq^q = 0. Recursively, we easily deduce that the approximate solutions exist 
if 

9.V = /, ^|y. =5, feC^m, g€C^{Y,) (5.53) 

has a solution. As ip{xi,X2) = f^^(^x2) ^ ~^ ^(^i) unique solution of (|5.53|) . this proves 

the existence of boundary layer expansions. These boundary layers are called Munk layers. 

5.2 Stommel and friction layers 

We now investigate the case where 

P ^ +00, > 0, Re, r given, (5.54) 

r 

which emphasizes the role of friction. In terms of e := (3^^ and := the last condition 
reads e^^^/ss 0. To fit exactly the context under consideration, we assume that Es = £^ for 
some rational number 6 < 1/3. The domain $7 remains the same as in previous subsection. 

Remind that the western symbol reads 
a«,(x',ei,6) := ie-H2 + 0(161 + 161) - ^e"' {fi + (1 + X^')6' 

+ 0(161 + 161)) - (e^ + i^ + x'^m + 0((|6I + 161)'; 

For all 6 £ the singular roots of a^(x',6) ■) satisfy either 

, /2W2 -1^1 '2\2t4 . . (-Re)V2 

j(l + X.)6 - + 6, ^■e■ 6 - [- 

Similarly, for the eastern symbol, 

-i I , (-Re)V2 
72 , or6 



There are two singular exponents, uniformly on Y: 71 = 5 {i.e. e'^'^ = Eg) and 72 = {i.e. 

The singular operators read 

a^,^i = iDg - rDg, a^.^^^ = -i Dg - rDg, 

1 2 1 2 

Oui,72 — ~fld — ——Dg, <^e,72 — —'''Id — ——Dg. 

ixe rte 
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Thanks to proposition El we deduce that e^'^ is a boundary layer size on l^.o and e'^^ is only 
a boundary layer size on Yf> (one root with positive imaginary part for a2?,e and , no such 
root for aj^)- 

We look for approximate solutions with 

nf„^ = if;^ + eT2-Ti^l + £2(72-71)^2 ^ 
Vbl = ^w,l ( 1 + V«;,2 ( ,^2 

, ,0 /^Xe(a;2) - \ 72-71,1 - x«,(a;2) \ , 
+ ^^6,2 ( j + ^^^^,1 I -^^ ,X2 1 + ... 

Equation yields (f^3^ . Equation yields 

^eV'^,2le=o = 5^, ^eV'e,2le=o = 9l, 

where 5^ g involve and "uj^ g 2 k < i — 1. All boundary layer operators have one root 
with positive imaginary part, thus only one condition can be given at the boundary. The last 
line prescribes Dgip^ e 2 ~ 0- Reasoning recursively, the existence of approximate solutions 
relies on the solvability of (|5.5:-{j) . The boundary layers with size Eg and ej := ^/ss/e are 
known as the Stommel and friction layers. 

5.3 Geostrophic degeneracy 

We consider again the asymptotics ()5.51|) . In subsection 15.11 we have considered a domain 
with western and eastern boundaries, but no northern or southern boundaries. Let us examine 
here the case of a closed basin. We assume simply that 

n = {xl + xl <l} , Y = dn = S\ 
Let us look for singular exponents near yo := (0,1). We chose local coordinates x'l = xi, 
1 — xf — X2- The local symbol reads, as |^i| + |^2| 00 



( It' 



V -6 / V -^-i"^ 



-^2 



Me?+(Y^ + i)e2^ + 0(161 + 161 

1 / 4x'2 



Let 6 £ Foi' y 7^ yoi i-G- x[ ^ 0, the roots 6 with |6| +00 satisfy 



which yields 7 = 1/3 in agreement with subsection 15.11 For y = i/q, i.e. x'^ = 0, this relation 
degenerates. The roots 6 going to infinity satisfy 

ze-i6 ~ ^^^2', ^■e.^2 ~ (i6Re)^/^e-i/4 
Re 
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so that 7 = 1/4. In particular, there are no singular exponents uniformly on a neighborhood 
of Ho = (0, 1) and (H3) fails. This phenomenon is known as geostrophic degeneracy. It takes 
place in various systems of fluid mechanics, including rotating fluids (as can be seen on figure 
n near the equator of the inner sphere, fohowing ^Hl)- To understand the structure of the 
solutions near such "turning points" is an open question, both from physical and mathematical 
viewpoints. 
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